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Abstract 

Solving the quantum-mechanical many-body problem requires scalable computational approaches, 
which are rooted in a good understanding of the physics of correlated electronic systems. Interacting 
electrons in a magnetic field display a huge variety of eigenstates with different internal structures, which 
have been probed experimentally in the Hall effect. The advent of high-performing graphics processing 
units has lead to a boost in computational speed in particular for classical systems. In the absence of a 
quantum-computer, it is thus of importance to see how quantum-mechanical problems can be cast into a 
seemingly classical dynamics, which can be efficiently implemented. At the same time, such mappings 
provide insights into the quantum-to-classical transition of many-body systems. 

1 Introduction 

The energy spectrum of interacting electrons in a magnetic field displays a rich structure, which results 
from the partial lifting of the degeneracies present in the non-interacting case. The magnetic field case bears 
some similarities with the ordinary harmonic oscillator, but the coupling of the angular momentum to the 
direction of the magnetic field leads to a reshuffling of the oscillator energies into highly degenerate Landau 
levels [1]. Once interactions between the particles are switched on, the oscillator states are no longer the 
eigenstates and in general numerical methods are required to obtain the energy-spectrum. Only for two- 
particles and specific Coulomb interaction values, exact analytic expressions are known [2 3 , 4|. The effect 
of interactions in a quasi two-dimensional gas of electrons at low temperature has received an enormous 
interest due to the observation of the quantized Hall effect, which shows up in measurements as plateaus of 
the transverse resistivity across a Hall bar for certain ranges of magnetic fields. The usual connotation of the 
integer and fractional Hall effects as being related to non-interacting and interaction electrons respectively, 
is somewhat arbitrary and questionable, since the electrons do not suddenly turn non-interacting in the 
integer Hall effect. Thus it is desirable to develop an understanding of both quantum Hall effects based 
on a common underlying Hamiltonian. Somewhat surprisingly, the classical Hall effect discovered in 
1879 [5 1 is presently not fully understood on a microscopic level. Recently, extensive numerical, first- 
principle calculations, became feasible using Graphics Computing Processors (GPUs), and have shown the 
importance of interactions and the proper incorporation of boundary conditions in the classical Hall effect 
in order to generate the observed Hall potential distribution [6]. These computations require large particle 
numbers (~ 10 4 ) to follow the self-consistent build-up of the Hall effect. The computed distribution of the 
Hall potential is in excellent agreement with measurements performed by Knott et al. Q and its possible 
implications for the quantum Hall effect are discussed in Ref. (SJ. In the converged solutions small regions 
next to the two contacts emerge, where the complete Hall potential drops over a tiny region in space. The 
continuum-solution based on conformal-mappings of the potentials contains a divergence of the electric 
field strength. The existence of such singularities points to non-trivial physics near the contacts, since a 
linear-response theory (based on zero bias-voltage) cannot spontaneously generate such points for a non- 
zero bias voltage. This motivates the search for possible similarities between the classical and quantum 
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Figure 1 : Part of the spectrum of three spin-less electrons in a magnetic field sorted by angular momentum. 
The quasi two-dimensional electrons interact via Coulombic forces. The labels denote mixed (m), sym- 
metric (s), and antisymmetric (a) symmetry of the spatial part of the wave function. For fermionic (spin 
1/2) electrons, the spatially symmetric states are not accessible. Energies are expressed in units of fvjji. 



Hall effects and highlights the need to understand the importance of interactions in all incarnations of the 
Hall effects. 

Here, for simplicity, I will focus on model systems with fixed electron numbers and not discuss transport 
properties. Already the computation of the quantum-mechanical few-electron case is a difficult many-body 
problem. Direct diagonalization methods in harmonic oscillator basis-sets are limited to small particle 
numbers, and interest exists to develop scalable methods from few to several hundred electrons. For the 
fractional quantum Hall effect, Laughlin developed a scalable method to sample information about specific 
many-body states. In the following, I review some known properties of Laughlin's wave-function and 
discuss how to use a classical map to compare the effect of the Pauli principle with the influence of the 
Coulomb interaction. 



2 Mapping of quantum mechanical systems to classical statistical 
mechanics 

One of my motivations is to discuss alternative methods to density functional theory (DFT), which employs 
a mapping of an interacting electron system to a non-interacting (still quantum-mechanical) Kohn-Sham 
system. A DFT implementation has to resort to approximate functionals to achieve this mapping, since 
the exact functional is unknown. In addition DFT is restricted to calculating ground-state properties and 
does not easily yield insight into the internal structure of the wave function. Recently, density-matrix 
renormalization group (DMRG) methods have been employed to analyze electrons in magnetic fields [9|. 
DMRG yields some information about excited states and the internal structure. A very different road is 
followed by methods which map the quantum-mechanical systems on classically looking equations. One 
example is coherent-state dynamics in combination with the time-dependent variational principle (TDVP), 
which allows one to derive systematically equations of motions for variational parameters of the wave 
function. The basic idea is to store the state of a many-body wave function in terms of trajectories of 
the variational parameters, which for coherent states include the positions, momenta, and possibly the 
widths of the time-evolving wave function J4][T0]. The stationary properties of the system are recovered 
by taking a discrete number of realizations along a trajectory and use these as a basis set for a subsequent 
diagonalization ifTTl . Laughlin proposed another approach, which maps a specific quantum state to a 
classical interacting system at finite temperature. In the following, I investigate how the incorporation of 
Coulomb interactions affects Laughlin's construction. 
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Figure 2: Pair correlation function of three electrons for the antisymmetric states shown in Fig.Q] from left 
to right panel L = 3,6,5 



3 The three electron case 

As first example, I discuss the spectrum of three electrons put in a plane perpendicular to a uniform mag- 
netic field. The spatial extension along the z-direction is neglected in the following and the Hamiltonian is 
given by 
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with Larmor frequency uii = #=■ and angular momentum L z ^ = —ih(xid V) i — yid Xj i), i = 1, 2, 3. The 
strength of the interaction is specified by k. For Coulomb interaction I set a = — 1, k > 0, and for repulsive 
oscillator interactions a = 2, k < 0. The Hamiltonian is separable into a center-of-mass part a relative part 
upon introducing Jacobi vectors 
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The center-of-mass part is contained in S3, and the relative part depends only on si, S2. 
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The introduction of Jacobi coordinates greatly facilitates the analysis of the three-body Hamiltonian. Espe- 
cially the analysis of possible symmetries can be linked to the Jacobi vectors, as has been demonstrated by 
P. Kramer and Moshinsky lfl2l and interaction matrix elements are simplified to integrals only involving 
two single-electron orbitals. Simonovic and Ciftja provide expressions for the exact Coulomb matrix ele- 
ments of few-body systems lfT3l[T4l . Laughlin also analyses the three-electron problem with Jacobi vectors 
in chapter 7 of Ref. [ 15 1. The diagonalization of the relative part results in the spectrum shown in Fig. Q] 
The angular-momentum of the relative motion is a good quantum number and for electrons in a magnetic 
field the distance between the electrons increases towards larger angular momenta. This in turn lowers the 
energy of the system. The number of states with a specific symmetry for each angular momentum can be 
derived from [ 12 1. lac hello and coworkers also analyzed algebraic approaches to the three-body problem 
in relation to molecular vibrations |[T6l[T7l . To gain insight into the internal structure of the states, I show 
in Fig. |2]the pair-distribution function 
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Figure 3: Lowest configurations of three electrons underlying the pair-distribution function shown in Fig. [2] 

for three antisymmetric states. The location of one electron is marked by the cross, while the shading 
indicates a high probability to find another electron. Two typical configurations underly these pictures, 
namely the triangular arrangement (Lagrange) and the colinear alignment (Euler) depicted in Fig. [3] 

4 Laughlin's plasma map 

Laughlin proposed that states for N electrons with coordinates Z{ = (xj + iyi)/lb of the form 

i> m (zi, . . . , z N ) = Y[(z 3 -z k ) m expl -^^2\ z i\ 2 \ ( 5 ) 

i<j I 1=1 ) 

are of special robustness against the effect of interactions. The magnetic length is given by l\, — -JH/ (2mcJi). 
The states are clearly antisymmetric under permutations of two particles and are states of angular momen- 
tum N ( N - 1 ) m _ These states are among the set of degenerate eigenstates of the non-interacting Hamiltonian 
(k = 0). Laughlin conjectured that the presence of interactions in general will (i) lift the degeneracies and 
(ii) will affect the family of states given by Eq. (0 less than all other states, thus leading effectively to a 
gap between the ground-state energy and the next excited state. The Jastrow factor depends only on the 
relative coordinates, while the sum in the exponent can be viewed as a center-of-mass term, which is sym- 
metric under exchange of two-particles. Laughlin illuminates the properties of these states by noting that 
the probability density can be written as a partition function: 

dz x ---dz N \^ m {z x ,...,z N )\ 2 = [ d Zl ---dz N e- 0u ^-< ZN \ (6) 
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j<k j 
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It is important to realize that this mapping to a classical system transforms the Fermi-Dirac statistics 
of the quantum-mechanical state at zero temperature to a classical Boltzmann system at temperature 
l/(fcsT) = 2m. The Pauli principle, which hinders two electrons to occupy the same state, is encoded 
in the logarithmic classical interaction of two point charges. A smooth density is achieved by translating 
the "fuzziness" of the quantum-mechanical location of the electrons to the thermal motion in the classi- 
cal ersatz-system. For instance, randomly distributed charges at initial time t = will start to move and 
collide and convert part of their electrostatic energy into kinetic energy. After a short while, an equilibra- 
tion between kinetic and electrostatic energies is reached, which allows one to determine the temperature 
of the system by looking at the histogram of the kinetic energy and to compare it to the Boltzmann dis- 
tribution exp(— (3Eki n ). Only for (3 = 2m do the classical density distribution n(r) and pair correlation 
function .(^(rl, r2) coincide with their quantum mechanical counterparts. Fig.|4]shows a visualization of 
the mapped state for m = 1, where electrons are at time t = randomly distributed and then a classical 
propagation using Newton's equations of motion is performed and the trajectories are plotted. The left 
panel shows the classical system at a too low temperature, while in the right panel the particles have the 
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Figure 4: Trajectories of the classical ersatz-dynamics for m = 1. Left panel: low temperature (T — 5 K), 
right panel: temperature T = 65 K corresponding to f3 = 2m. By evaluating the histogram of the particle 
positions and distances, the density and pair correlation function of the equivalent quantum-mechanical 
state are evaluated. 



E = 12.7 X 1 








E = 12.7 12 






c 


t 






> - 








s 

























Figure 5: Pair correlation function for three Coulombic -interacting electrons in a magnetic field. Left 
two panels: weak Coulomb interaction, left-most panel: exact diagonalization, second panel histogram of 
the pair-distribution of the classical ersatz-dynamics using Eq. ©. Right two panels: increased Coulomb 
interaction, exact diagonalization vs. ersatz-dynamics (right-most panel). 



correct temperature to reproduce the statistics of |-0i (zi, . . . , Zigg)| 2 . Several methods are available to de- 
termine the pair-correlation function: Laughlin used the hypernetted-chain method [ 15] and the Metropolis 
algorithm is discussed in Ref. ifTSl . Here, I will use a direct integration of the classical equations of mo- 
tions and obtain the density and pair correlation function by sampling the electron positions periodically. In 
the classical ersatz-dynamics, the effect of the Coulomb interaction is included by amending the classical 
potential containing the logarithmic "Pauli-potential" with a Coulombic term: 




In order to fill out phase-space for the three electron case, an ensemble average over different initial posi- 
tions of the electrons is performed. In Fig.|5]l compare the pair correlation function of the mapped system 
with the result from an exact diagonalization of the Hamiltonian. The main effect of the Coulomb interac- 
tion is a stretching of the pair-correlation function, seen as an increase in the relative distances between the 
electrons. 
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Figure 6: Radial density profile of an 196 electron quantum dot, including a disk-shaped positive back- 
ground charge. The dashed line shows the non-interacting case, which represents the lowest Landau level 
(to = 1). The solid line denotes the result including the Coulomb term in Eq. |9). DFT results for similar 
systems are shown in |fl9l . 

5 Larger systems 

The three-electrons in a magnetic field were an important test-case for Laughlin's development of a theory 
of the fractional quantum Hall effect. It is desirable to investigate how the classical ersatz-dynamics works 
for larger particle numbers. This has been investigated by Monte-Carlo samplings fl8l . In addition to the 
Coulomb forces, in this case also a uniform, disk-shaped, positive background is taken into account. Here, 
I present results using again a direct integration of Newtons equations for 196 electrons, which are in good 
agreement with the results shown in Ref. IfTSIl . The integration is performed very efficiently using GPU 
computers. Similar systems have been investigated using current density-functional theory |fT9l l20ll . 

6 Summary 

I have reviewed some aspects of electrons in a magnetic field. The classical ersatz-dynamics provides an 
intriguing representation of the Pauli principle in the form of an interacting electron gas at finite tempera- 
ture. The addition of Coulomb forces yields good results for Laughlin's states already for three electrons. 
For more general electronic systems, the classical-ersatz dynamics has seen some renewed attention, see 
for instance the overview in Ref. [21|. If the ersatz-dynamics provides a viable alternative to other many- 
body methods such as DFT, DMRG, and TDVP, requires further investigation and a careful comparison of 
the results for different systems. One difficulty of Laughlin's map is its reliance on a specific state, rather 
than the ability to obtain a whole set of eigenstates. From a computational point of view, the classical dy- 
namics can be excellently implemented on the emergent super-computing hardware consisting of graphics 
processing units. 
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